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Abstract. In the paper, we consider the large time behavior of solutions to the 
convection-diffusion equation ut — An + V • f{u) = in iR" x [0, oo), where /(u) ~ u'^ 
as u — > 0. Under the assumption that q > 1 + l/{n + jS) and the initial condition 
uq satisfies: uq G L^{]R^), Jjp^„ uo{x) dx = 0, and ||e*^uo||Li(jR") < Ct~^/'^ for fixed 
P G (0, 1), all t > 0, and a constant C, we show that the L^-norm of the solution to 
the convection-diffusion equation decays with the rate as t — > oo. Moreover, 

we prove that, for small initial conditions, the exponent q* = l + l/{n + P) is critical 
in the following sense. For q > q* the large time behavior in LP{]R^), 1 < p < oo, 
of solutions is described by self-similar solutions to the linear heat equation. For 
q = q*, we prove that the convection-diffusion equation with f{u) = u\u\'^*~^ has 
a family of self-similar solutions which play an important role in the large time 
asymptotics of general solutions. 
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1 Introduction 

In this paper, we study the large time behavior of solutions u = u{x, t) [x G 
iR*^, t > 0) to the Cauchy problem for the nonlinear convection-diffusion 
equation 

(1.1) ut- An + a- V(m|m|'?-1) = 0, 

(1.2) m(x,0) = no(x), 

where g > 1 and the vector a G iR" are fixed. The assumptions Uq G L^{IR^) 
and /jfjn uq{x) dx = Q will also be required. 

The typical nonlinear term occurring in hydrodynamics in the one-dimension- 
al case has the form uu^ = {v? /2)x (as in the case of the viscous Burgers 
equation). The most obvious generalization of this nonlinearity consists in re- 
placing the square by a power where g is a positive integer. Here, however, 
we intend to observe a more subtle interaction of the nonlinearity with dis- 
sipation, consequently, we need to consider a continuous range of parameters 
q. The problem then appears with the definition of u'^ for negative u and for 
non-integer q. In order to avoid this difficulty, we chose the nonlinear term of 
the from a ■ 'W{u\u\^~^). This was done only to shorten notation in this report. 
Note that, in fact, the following property of the nonlinearity will be essential 
throughout this work: 

• the nonlinear term in ( |1.1| ) has the form V ■ f{u) where the C^-vector 
function / satisfies \ f{u)\ < C|m|'', |/'(m)| < C|m|^~^ for every u & M, 
q > 1, and a constant C. Moreover, if the balanced case is considered 
(i.e. g = 1 + l/(n + /5)), the limits 

lim f{u)/\u\'^, and lim f{u)/\u\'' 

should exist and the both should be different from 0. 

Recent publications developed versatile functional analytic tools to study the 
long time behavior of solutions of this initial value problem. 
Concerning the decay of solutions of (|l.l| )-( pr2|) and, more generally, of scalar 
parabolic conservation laws of the form Ut — Au + V ■ f{u) = with integrable 
initial conditions, Schonbek pOf was the first who proved that the L^-norm 
tends to as t — s> oo with the rate To deal with this problem, she 

introduced the so-called Fourier splitting method. The results from [E9| were 



extended in the later work |3^, where the decay of solutions in L^(iR"), (1 < 
p < oo) was obtained, again, by a method based on the Fourier splitting 
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technique. It was emphasized in |30| that the decay rates are the same as for 
the underlying hnear equations. 

Next, Escobedo and Zuazua proved decay estimates of the L^-norms of so- 
lutions by a different method under more general assumptions on nonlinearity 
and under less restrictive assumptions on initial data. Finally, by the use of 
the logarithmic Sobolev inequality, Carlen and Loss ^ showed that solutions 
of viscous conservation laws satisfy 



^^0 



for each l<r<p<cxo, allt>0, and a numerical constant C > depend- 
ing on p and q, only. Here, we would also like to recall results on algebraic 
decay rates of solution to systems of parabolic conservation laws, obtained by 
Kawashima [^, Hopf and Zumbrun [^, Jeffrey and Zhao and Schonbek 
and Siili [^. Smallness assumptions on initial conditions were often imposed 
in those papers. 

The first term of the asymptotic expansion was studied as the next step in 
analysis of the long time behavior of solutions to ( |1.1|) - (|1.2|) . Assuming that 
Uq e L^{1R'^), roughly speaking, these results, cf. e.g. [|, |14|, [1^, [11], 

I, a I, P, |2l, fall into three cases: 



Case I: q > 1 + when the asymptotics is linear, i.e. 

(1.3) _ MG{-,t)\\p ^ as t ^ oo, 

where M = J]i^r^uo{x) dx, G{x,t) = (47rt)""/^ exp(— |xp/(4t)) is the fun- 
damental solution of the heat equation. Hence, this case can be classified 
as weakly nonlinear, since in this situation the linear diffusion prevails 
and the nonlinearity is asymptotically negligible. 

Case II: q = 1 + when 

(1.4) t("/2){i-i/p)||^(.^^) _ UM{-,t)\\p ^ as t oo, 

where UM{x,t) = t~"'/'^UM{xt~^^'^ ,1) is the self-similar solution of ( [1.1|) 
with Uo{x) = M6o. Here, diffusion and the convection are balanced, 
and the asymptotics is determined by a special solution of a nonlinear 
equation. 

Case HI: l<g<l + l/n, when the convection points in the x„-direction 
(i.e. a= (0, ...,0, 1)). Here 



(1.5) _ UM{-,t)\\p ^ as t 



oo, 
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holds, where Um is a particular self-similar solution of the partly viscous 
conservation law Ut - Ayll + ^{U\U\'^-^) = such that Uo{x) = MSq 
in the sense of measures. Here x = (y, y = (xi, . . . , and 

= I^i=i Hence, the asymptotics of solutions is determined by 

j 

solutions of an equation with strong convection and partial dissipation. 



Finally, we recall that, in the weakly nonlinear case, Zuazua |^ found, for 
solutions to (|1.1|) - (|1.2|) , the second order term in the asymptotic expansion as 
t — > cxo. He observed that asymptotic behavior of the solution differs depending 
if q satisfies 1+1/n < q < 1+2/n, q = 1+2/n, or g > 1+2/n. Analogous results 
for Levy conservation laws were obtained in @, ^, and for convection-diffusion 
equations with dispersive effects in p2|, |2^. Related results on the stability 



in L^(iR") of traveling waves (or shock waves) in scalar viscous conservation 



laws can be found in the papers by Serre [^] and Freistiihler and Serre |]T7 
Some results on the L^-stability of the zero solution of degenerate convection- 
diffusion equations can be found in the article by Feireisl and Laurengot |T^ . 



Here, we assume that M = Jjiin UQ^x) dx = Jji^nu{x,t) dx = 0, thus the 
corresponding self-similar intermediate asymptotics in (|1.3|) - (|1.5| ) are equal to 
for every g > 1. Moreover, for p = 1 the asymptotic formulae in ( |1.3|) -( pT5D 
say nothing else but ||'u(-, — > as t — ^ oo. 

The goal of this paper is to find self-similar asymptotics in L^{1R") of solutions 
to ( p..lD -( |]r2|) with M = imposing additional conditions on the initial data. 



We assume that uq satisfies ||e*^Mo||i < Ct'^^"^ for some P G (0, 1), all t > 0, 
and C independent of t. Such a decay estimate of solutions to the linear heat 
equation is optimal for a large class of initial conditions (cf. Propositions 
|2J] and |3.1| below). Under these assumptions, we improve the known algebraic 



decay rates of the solutions to (|1.1|)-(1.2) in the L^-norms for every 1 < p < oo. 



In addition, if the initial data are sufficiently small, we discover the new critical 
exponent g* = 1 + l/(n + /5) such that 



for q > q* the asymptotics of solutions to (|1.1| )- (|1.2| ) is linear and de- 
scribed by self-similar solutions to the heat equation (cf. Corollaries |2]2 
and |3]5, below) ; 

q = q* corresponds to the balanced case, and the asymptotics of solu- 
tions corresponding to suitable small initial conditions is described by 
a new class of self-similar solutions to the nonlinear equation ( p..l| ) (cf. 



Theorem 2.3 and the discussion in Section 5. 
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In the next section of this paper, we briefly present the main results. The 
results corresponding to the case q > q* are contained in Section 3. Section 4 
considers the case when the exponent q = q* is critical. In Section 5, we explain 
how to derive, from our general theorems, self-similar solutions to the nonlinear 
equation ( |1 . 1|) and how to study large time asymptotics of general solutions. In 
the last section, we discuss possible applications of our ideas to other equations 
such as the Navier-Stokes equations, the KdV-Burgers equation, and the BBM- 
Burgers equation. 

Notation. The notation to be used is mostly standard. For 1 < p < oo, 
the L^-norm of a Lebesgue measurable real-valued function defined on iR" is 
denoted by ||f We will always denote by || ■ Ha" the norm of any other Banach 
space X used in this paper. 

If /c is a nonnegative integer, W^'''^{1R") will be the Sobolev space consisting 
of functions in L'p^IRJ^) whose generalized derivatives up to order k belong to 
U'{]BJ'). 

The Fourier transform of v is defined as v{^) = (27r)^"/^/^„ e~^^^v{x) dx. 
Given a multi- index 7 = (71,..., 7„), we denote d'^ = d'^'^^ / d1\...d1"- On the 
other hand, for /? > 0, the operator is defined via the Fourier transform as 

(i5%;)(o = mm- 

The letter C will denote generic positive constants, which do not depend on t 
and may vary from line to line during computations. 



2 Main results and comments 



We recall that for every uq G L^(iR"), the Cauchy problem ( |1.1|) - (|1.2|) has a 
unique solution in C([0, cxo); L^(iR")) satisfying 

u e C((o, 00); iy2,P(iR")) n c^((o, 00), (iR")) 

for all p G (1,00). The proof is based on a standard iteration procedure 
involving the integral representation of solutions of ( |l.l| )-( p^ 



(2.1) 



u{t) 



a- Ve(*~")^(M|Mr-^)(r) dr 



19-1^ 



(see, e.g. [|T3| for details). Here, c^^Uq is the solution to the linear heat equation 
given by the convolution of the initial datum uq with the Gauss- Weierstrass 
kernel G{x,t) = (47rt)~"/^ exp(-|xp/(4t)). Formula (|2l]) will be one of the 



main tools used in the analysis of the long time behavior of solutions. 
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Let us also recall that sufficiently regular solutions of (1.1)-( [L^) satisfy the 
estimate 

(2.2) h(-,t)||,<C(p,r)r("/2)(i/'-Vrt||^^||^ 

for all 1 < r < p < oo, all t > 0, and a constant C{p,r) depending on p and 
r, only. Inequalities (|2.2|) are due to Carlen and Loss H, Theorem 1]. We 
also refer the reader to |^ where counterparts of ( |2.2|) were proved for more 
general equations: so-called Levy conservation laws. 

Section 3 contains the analysis of the large time asymptotics of solutions to 
the linear heat equation. Easy calculations show that for every mq G L}{]R"') 
such that /jRuMo(a;) (ix = we have ||e*^Mo||i — as t ^ oo. The following 
proposition asserts the existence of a large class of initial conditions for which 
the large time behavior of e*^Mo is self-similar. Here, we need the notion of 
the Riesz potential Ip and the fractional derivative defined in the Fourier 
variables as 

(2.3) (V^)(0 = ^ and (i5^t.)(0 = I^I^^O- 

Proposition 2.1 Let /5 > and 7 = (71, 7„) he a multi-index with 7i > 0. 
Assume that I^Uq G L}{IR^). Denote 

(2.4) ^ = ,lrnW= / (Vo)(a:)rfa;. 
Then 

(2.5) ||9V^wo||i<C^t~^/'~l"l/l//?wo||i 

for allt > and C = C{f3,'y) independent oft and uq; moreover, 

(2.6) t^/2+|7l/2||^7etA^g(.) _ Ad^D^G{-,t)\\i 
as t ^ 00. 



Let us emphasize that we do not assume that /5 < 1 in Proposition |2.1|. This 



condition only becomes necessary when the convection term is present. Here, 



we also refer the reader to Proposition where self-similar asymptotics of 
the heat semigroup is studied under more general assumptions on initial data. 

In our first theorem on the large time behavior of solutions to the nonlinear 



problem (|l.l|)-(pr2|), we assume the decay of ||e MoHi with a given rate and we 



prove that the same decay estimate holds true for solutions to (|Ll|)-(pr!^). 
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Theorem 2.1 Fix < (3 < 1. Assume that uq E L^{M') (1 L'i{IBJ') satisfies 
the inequality 

(2.7) ||e*^Mo||i < Cr'^/' 

for allt > and a constant C independent oft. Let u he the solution to (U^Jj- 
l\1.2j with Uq as the initial datum. If q> 1 + 1/n, then there exists a constant 
C such that 

(2.8) \H;t)\\i<C{l+t)-^/' 

for allt > 0. The estimate p.<S| j holds also true for 1 + 1/ (n+P) < q < 1 + 1/n, 
with < f3 < 1, provided uq G ^^(iR") n L'=^{IRJ') and sup^^o ^^^^l|e*^wo||i «s 
sufficiently small. 



Remark 2.1. The assumption ( p.7|) means that uq belongs to the homoge- 
neous Besov space Bi^'^ (cf. ( |2.11| ), below) which will play an important role 
in the analysis of the balanced case q = 1 + 1/ {n + (3). □ 



The approach formulated in Theorem ^T!], saying that the decay estimates 
imposed on the heat semigroup lead to the analogous estimates of solutions to 
a nonlinear problem, appears in several recent papers. Here, we would like only 
to recall (the list is by no mean exhaustive) the works on the Navier-Stokes 
system by Schonbek [RTl and Wiegner p6[ where the L^-decay of solutions 



was studied as well as by Miyakawa PB[ where decay of the L^-norm and IHP 



norms (the Hardy spaces) of weak solutions was shown. Moreover, our results 
extend essentially the recent paper by Schonbek and Siili where general 
conservation laws were considered. 

If we combine the decay from ( ^.8] ) with inequalities ( [^.2D , we obtain the im- 



proved L^-decay of solutions to ( |l.l| )-(p^). Moreover, applying such estimates 
to ( |2.1|) we find the asymptotics of solutions for q > 1 + l/{n + (3). The 
following corollary contains these results. 



Corollary 2.1 Under the assumptions of Theorem \2. 1\ , for every p G [1, oo] 
and (3 G (0, 1), there exists C = C{uo,p) independent oft such that 

(2.9) ||M(-,t)||p < C(l + t)-("/2)(^-i/'')-^/' 

for all t > 0. Moreover, for q > 1 + l/{n + (3) and for every p G [1, oo] it 
follows 

(2.10) t("/2){l-lM+/3/2||^(.^^)_giA^^(.)||^^0 t^oo. 
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A slightly stronger version of this corollary is formulated and proved in the next 
section (cf. Corollary |3.1| , below). Here, we only emphasize that combining 
(p.lO|) with Proposition |2.1| we obtain that the large time behavior of solutions 
to (|1.1|) - (|1.2| ) with q > 1 + 1/n (or, if the data are sufficiently small, for 
q > 1 + l/{n + /3)) is described by special self-similar solutions to the heat 
equation. This is worth stating more precisely. 



Corollary 2.2 Under the assumptions of Theorem \2.1\ and Proposition 



(or Proposition |g. j| with = |^|^, see Section 3) the solution to 
with q > 1 + l/{n + (3) satisfies 

tHm-m+mi\u^.^t)-AD^G{;t)\\p^O as t^oo. 

Our next results, studied in Section 4, correspond to the balanced case 

1 



q = q 



n + p 

for some fixed < (3 < 1. We will work in the homogeneous Besov space 

B^^^°° = {ve 5'(iR") : ||t;||g-/3,oc < oo}, 
where S'{1R") is the space of tempered distributions and the norm is given by 
(2.11) M.-P,oo = sup/Z^e'^t^lli. 

1 s>0 

The standard way of defining norms in Besov spaces is based on the Paley- 
Littlewood dyadic decomposition. The choice of the equivalent norm ( p. 11 



allows us to simplify several calculations. Recall here that Proposition |2.1| 
describes a large subset in of initial conditions Uq. 

Section 4 contains the proofs of two main theorems. The first one provides a 
construction of global-in-time solutions to (|l.l|)-(pr^) with q = l + l/(n + /?) 



and suitably small initial data in the space i3f The second theorem gives 
asymptotic stability of solutions in the balanced case. The precise statement 
of the theorems is the following. 

Theorem 2.2 Fix (3 G (0, 1) and put q = 1 + l/{n + (3). There is e > Q such 
that for each uq G satisfying H-uoHg-^i.oo < e there exists a solution of 

(f7j|j-(i77|j for allt>0 in the space 

X = C([0, oo) : i3r^'°") 

n {u : (0,cx)) ^ L«(iR'^) : sup t("/2){i-i/g)+/3/2 ||^ < 



Convection- diffusion equation 

This is the unique solution satisfying the condition 

,^W2)(l-l/,)+/3/2|, .^^|| <2£. 



sup V 1 1 o y 1 1 



Theorem 2.3 Let the assumptions from Theorem ^.Sj hold true. Assume that 
u and V are two solutions of j| j- pT^ constructed in Theorem corre- 
sponding to the initial data uo,vq G respectively. Suppose that 

(2.12) \imt^/'\\e'^iuo-vo)\\i = 0. 

t — >CXD 



Choosing e > in Theorem |^.^ sufficiently small, we have 



(2.13) lim t("/2)(i-i/p)+/3/2||^(.^^) _ = 



for every p E [1, oo]. 



In Section 5, we show how to use Theorem p.2| in order to obtain self-similar 
solutions to equation ( |1.1| ) with the critical exponent q = l + l/{n-\- (3). More- 
over, we explain the role of self-similar solutions in the large time asymptotics 
of other solutions to (|l.l|)-(0 



3 Asymptotics of solutions for q > l + + 

As noted in Section 2, the first problem is to find a class of data that will 
insure the decay of solutions to the heat equation in L^(iR"). This is obtained 
in Proposition |2.1j , where this class of data is shown to be constituted by 
functions such that their convolutions with Riesz potentials lie in L^(iR"). 
Now, we establish Proposition 

Proof of Proposition |2.1J . Let us note that the limit in (|2.4| ) exists, since 
^o(0/l^l^ is continuous as the Fourier transform of an integrable function I/^uq. 
First, we prove that 9^L'^G(-, 1) G L^fW"). Obviously, d^Dt^G{-, 1) is bounded 
and continuous because its Fourier transform (^O'^^l^l^^"'^'^ is integrable. More- 



over, it follows from Ch. 5, Lemma 2] that for every (3 > there exists a 



finite measure fip on given by 



Mi) 
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Hence, d'^D^G{-, 1) = /i/3 * where the function Kp^^ is defined via the 
Fourier transform as Kp^^{^) = + \^\'^)^/'^e-\^\\ It is easy to prove that 

Kp^^ G iS(iR") (the Schwartz class of rapidly decreasing smooth function), and 
this implies the integrabilty of d'^D^G{-, 1) for every multi-index 7. 
Now, the change of variables yields that d''D^G{x, t) has the self-similar form: 



d^D^G{x,t) = t-"/2-/3/2-|7|/2^^7/)/3G')(2;/v^, 1) 



(3.1) 



for all X e ]RJ^ and t > 0. 

To prove ( p.5|) , use the Young inequality for the convolution, and thus by ( p.l|) 
it follows 

Wd^e'^uoWi = \\d^D^G{t)*If3Uo\\i 

< ||9^D^G'(-,t)||i||Vo||i 

< r^/'-i^i/l9^/^''G(-,i)||i||Vol|i 

for all t > 0. 

For the proof of (|2.6|) , observe that the change of variables z = x/ \/t combined 
with ( p. ID leads to the following expression 



dx 



(3.2) =t/^/2+H/2 



d^D(^G{x -y,t)- d^D/^Gix, t)] IpUo{y) dy 



< 



\IpUoiy)\ id^'Df'G){z-y/Vt,l)-{dWf'G){z,l) dydz 



From the first part of this proof, the function d^D^G{z, 1) is continuous, hence 
the integrand on the right hand side of (|3.2|) tends to as t — oo for all 
y, z G iR". Denote 

A{z,y,t) ^ {d^DPG){z - y/Vi, 1) - {d^D^G){z, 1). 

To apply the Lebesgue Dominated Convergence Theorem to the integral on the 
right hand side of ( p.2| ), it is necessary to show that there exists F G L}{]R"') 
independent of ?/ G iR" and t > 1, such that 

(3.3) \A{z,y,t)\<F{z) 

for all 2, y G iR" and t > 1. Note that 



A{z,y,t) 



-iy/Vi _ 
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Moreover, the symbol b{^,y,t) = {1 + \p-iyl^i - i] e-l^l' is a C= 

function of (^,1/) G iR" x iR", and satisfies the differential inequalities 



\dldlh{i,y,t)\<C{a.l.^)iX^\i\) 



-N-a 



for all multi-indices a and 7, all G W, and C{a,'y,N) independent of 
E and t > 1. By ||3^, Ch. VI, Sec. 4, Prop. 1], the (inverse) Fourier 
transform with respect to ^ of b{^, y, t) satisfies the estimate 



\J^^%-,y,t){z)\<C{l + \z\)- 



■N 



for all G -ffV, and a constant C = C{N) independent of z,y & IR^ and t > 1. 
Finally, the use of the measure /^^ from the first part of this proof combined 
with standard properties of the Fourier transform and the convolution lead to 
the representation A{-, y, t) = fxp* T^^h{-, y, t). Hence, ( p.3| ) holds true for the 
function F{z) = C[fip * (1 + | ■ |~^)](z) with any N > n. This completes the 
proof of Proposition |2.1|. □ 



We recall that, in Miyakawa obtained the L^-decay of e*^Mo provided 
the |xp-momentum of the data is bounded. Below, we will show that our 
assumptions is weaker than the one assumed by Miyakawa. 

Remark 3.2. The L^-decay of solutions to the linear heat equation formulated 
in ( ^.51 ) was proved by Miyakawa under the assumptions 

(3.4) Mo G L^(iR"), / uoix)dx = 0, / \x\'^\uoix)\ dx < 00 

for some < f3 < 1. To show that our assumption //3M0 G L^{1R'^) is weaker 
than (|3.4]) it suffices to establish the inequality 



(3.5) llVolli < C / \xf\uo{x)\ dx 

JlR" 



valid for every uq satisfying ( |3^ ) with /? G (0,1). Let us sketch the proof 
of (P3|), however, it does not play any role in our considerations, below. It 
is well known that (//3Mo)(x) = C(/3,n) /jfj„ \x — y\'^~"'Uo{y) dy (in fact, this 
representation holds true for every /3 G (0,n)). Hence, using the assumption 
/jR" '^oiy) dy = and changing the order of integration yield 



\\If,uo\\i<CiP,n) 



dx \uo{y) \ dy. 



Next, note that the integral with respect to x in the inequality above is finite for 
every y G IR^, because its integrand ||x — — is locally integrable 
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and behaves like Ix]'^ ""^ " as |x| oo (here, the assumption /5 G (0, 1) is 
crucial). Hence, by the change of variables, it follows that 



1 



\x-y 



\X 



\n-f3 



dx 



\y\ 



iR" 



1 



\uj -y/\y\\'"'-(^ 



duj. 



Since supygjjn^jo} Jmn \ \u! — y/\y\\'^ " — \uj\^ "| du < oo (the proof of this ele- 
mentary fact is omitted), we obtain ( p.5| ). □ 

The self-similar asymptotics of e^^uo in L^^M"') with p G [2, oo] can be derived 
under weaker assumptions on uq. This is stated in the following proposition. 

Proposition 3.1 Leti = i{C,) denote a function homogeneous of degree (3 > 0. 
Assume that uq satisfies 



(3.6) 



«o(0 ^ 
CeK"\{0} «(4j 



and 



^o(0 

liPo £(0 



lim 



A 



for some A & M. Denote by C the Fourier multiplier operator defined via the 
formula Cv{^) = i{^)v{C,). Under these assumptions, for every p G [2,oo] and 
for every multi-index 7, it follows 

^n(l-l/p)/2+/3/2+H/2||^7etA^^_^57£G'(t)||p^0 aS t ^ OO. 

Proof. The main tool here is the Hausdorff-Young inequality 
(3.7) \\v\\p<C\\v\\,, 



valid for every l<q<2<p<oo such that l/p+ 1/q = 1. Hence, (|3.71 ), the 
change of variables ^t^^'^ = co, and the homogeneity of £ yield 



Id^e'^^uo- Ad''CG{t)\\l 



< C 



p 

7p-t|5l = 



( u,{i)-Am \ 

\ m ) 



^^_n/2-(/3/2+|7|/2)g 



di 

e{uj/ty^) 



A 



dC 



Now, the assumptions on uo in ( |3.6| ) allow us to apply the Lebesgue Dominated 
Convergence Theorem in order to prove that the integral on the right hand 
side tends to as t — cxd. □ 



Remark 3.3. The conditions formulated in ( p.6| ) appear in a natural way 
if Hardy spaces are considered. Let us recall that a tempered distribution 
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V belongs to the Hardy space on iR" for some < p < oo whenever 
v+ = supi>o \{(pt*v)\ G LP{1R^), where (f)t{x) = t~''(f){x/t) with G 5(iR") 
such that /jjn dx = 1. We refer the reader to |3^ where several prop- 
erties of Hardy spaces are derived. We recall that is a Banach space 
strictly contained in L}{]R^) and that U'{]Fr') = Ti^ for p > 1 with equiv- 
alent norms. Suppose now that p < 1 and Uq G 'HP. It is known (cf. |3^, 



Chapter HI, §5.4]) that the Fourier transform uq is continuous on iR" and 
|mo(OI ^ C'|^|"(^/P~^)||uo||-HP for all ^ G iR". Moreover, near the origin, this 
can be refined to \im^^QUQ{^)\^\~'^^^/P~^'> = 0. Hence, assumptions ( |3.6| ) are 
satisfied with = |^|^, p G (0, 1), and A = 0, if e.g. Uq G 7t:"/("+/3). □ 



Theorem jO] is the main decay theorem proved in this section. It ensures that 
the L^-norm of solutions to the convection-diffusion equation decay at the rate 
^-/3/2 provided their initial data are such that the corresponding solutions to 
the heat equation decay at the same rate. 



Proof of Theorem |2.1| . The proof of this theorem relays on a systematic 



combination of the integral equation ( |2.1D with inequality ( |2.2| ). Note that 
since Uq G ^^(iR") n L«(iR"), by (^), it follows that 

(3.8) \H;t)\\l < C(ho||i, lko||,)(l + t)-("/2)(^-^) 

for all t > 0. Hence, computing the L^-norm of (|2.1|) , using the assumption on 
Mo, and (|3.8|) yield 



|^^(-,t)||i < \\e''^uoh+ f'\\a-VGi-,t-r)\\i\\u{;t)\\ldT 



q 

t 



(3.9) < + c f {t- r)-i/2(i + ^)-W2)(.-i) dr 

Jo 

< Ct~^/^ + C{ t~'/Hog{e + t), for g = 1% |; " 

rV2, for g > 1 + i 



For g > 1 + (l3 + l)/n, estimate (|2.8|) follows immediately from ( p.9|) , since 
-1/2 < -(3/2 and 1/2 - (n/2)(g - 1) < -/5/2 for this range of g. 
Next, consider l-|-l/'ri<g<l-|-(/3-|-l)/n. A simple calculation shows that 
a = —(1/2 — {n/2){q — 1)) satisfies < a < (3/2. Moreover, it follows from 
(^ that 

(3.10) ||M(-,t)||i < C(l + t)"". 
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Combining inequality (|2.2| ) with ( 3.101) yields the improved decay of the L'^- 
norm 

(3.11) ||^"(-,t)||, < C(l+t/2)-("/2)(i-i/.)||^(.^^/2)||^ 

Hence, repeating the calculations from ( p.9|) , using ( p.ll|) instead of ( p.8|) , gives 

(3.12) ||M(-,t)||i < At-r)-i/2(l + r)-("/2)(^-i)-'?"rfr. 



If — (n/2)(g — 1) — ga < —1, the integral on the right hand side of (p.l2|) tends 
to as t ^ oo faster than t~l^^'^ and this ends the proof. On the other hand, 
if —{n/2){q — 1) — qa > —1, by the definition of a, it follows from ( p.l2| ) that 

||M(-,t)||i <cr^/' + cr°(^+^). 

Hence, if —a{q + 1) < —f3/2, the proof is complete. If, on the contrary 
—a{q + 1) > —(3/2, we have the new estimate 

lk(-,t)lli<c(i + t)-°(«+^), 

which we use as in ( |3.11| ) (with a replaced by a(g + 1)) to get an improved 
decay of the L^-norm: \\u{-,t)\\g < C{1 + t)-W2)(i-i/g)-"(g+i). Consequently, 
a finite number of repetitions of the above steps yields (p.8|). 
Finally, let us prove ( p.8[ ) for l + l/{n + l3) < q < 1 + 1/ n under the assumption 
that supj^g ||e*^Mo||i is sufficiently small. For simplicity of notation, we put 

g = g* = i + ^-^, 

n + P 

and we use systematically the following inequality (obtained from the Holder 
inequality and from (p.2|)) 



(3.13) \\ui.,t)\\i < im-,t)ii^:iK,t)ii^-^* < c{\\uo\u\H.,t)\\i: 

for all t > 0. To proceed, we also define the auxiliary nonnegative continuous 
function 

g{t)^ sup p/2||^(.,r)||i)+ sup M^/'+^Z^)/''* h(-, r)l|,,) . 

Now, computing the L^-norm of the integral equation ( p.l| ) and using ( |3.13|) 
yield 

t^/'M-Mi < t'^/1e*^.o||l + C^^/^/J(^-r)-Vl..(,r)||^ dr 

(3.14) < t^/^e'^Molli 
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for all t > 0. An elementary calculation shows that the quantity 





is finite for every t > (since < (3 < 1) and independent of t. A similar 
reasoning gives 

(3.15) \\u{;t)\\,, < (t/2)-("/2)(i-i/'?*)||e(*/2)A^^||^ 

+ g'^'it) C [\t - r)-W2)(l-V5*)-l/2^-l/2-/3/2 

Jo 

Note that -{n/2){l - 1/g*) -/5/2 = -(1/2 + /3/2)/g*. Moreover, the quantity 

^(l/2+/3/2)/q* f^^^ _ ^^|-(n/2)(l-l/g*)-l/2^-l/2-/3/2 







is finite (since — (ri/2)(l — 1/q*) — 1/2 > —1) and independent of t (by the 
change of variables). 

Combining inequalities ( |3.14|) and (|3.15| ) yields 



(3.16) g{t) < Cisup Tf^/^\\e'^uo\\i + 029''* (t) 

0<T 

for all t > and constants Ci and C2 independent of t. 
Finally, let 

F{y) = A + CiV"' - y where A = Ci sup T^^'^\\e'^uo\\i 

0<T 

and where g* > 1. If A > is sufficiently small, there exists yo > such 
that -F(|/o) = and F{y) > if ?/ G [0,yo). Moreover, it follows from ( p.l6| ) 
that F{g{t)) > 0. Since g{t) is a nonnegative, continuous function such that 
g{0) = 0, we deduce that g{t) G [0,yo) for all t>0. This completes the proof 
of Theorem 2.1. □ 



As a consequence of Theorem |2.1| we get Corollary |2.1| . Actually, here we prove 
its slightly stronger version. 

Corollary 3.1 Under the assumptions of Theorem for every p G [1, cxd] 
and (3 G (0, 1) there exists C = C{uq,p) independent oft such that 



(3.17) 



M(-,t)||p<C(l + t)-("/2){i-iM-/3/2 
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for all t > 0, and 

||n(-,t)-e*^«o(-)llp 



(3.18) 



r("/2)(i-Vrt-i/2 iog(e + t) for q = ^, 

^-{n/m-l/p)-l/2 J g > % 



for allt>l. 



Proof. Inequality ( |3.17|) is obtained combining ( p.2|) with ( p.8|) as in ( p.ll|) 
where q is replaced by p and a by /?/2. 



In view of the integral equation (^.11), to prove ( ^.18| ) it suffices to estimate 



the L^-norm of the second term on the right hand side of ( [2.11 ). Here, split 
the integration range with respect to r into [0,t/2] U [t/2,t] and study each 
term separately as follows. Using the Young inequality for the convolution and 
( CT) yields 



ft/ z 

/ ||a- Ve(*-")^(u|u|'?-i)(r)|Lrfr 
Jo 

(3.19) < J^^'\\a-VG{;t-T)Uu{;T)\\ldT 



rt/2 

< C _ ^)-W2)(l-l/p)-l/2(^^ _^ ^^)-(n/2){g-l)-/3<^/2 



/O 



< C{ r("/2)(i-Vp)-i/2 iog(e + for g = ^, 

^-(n/2)(l-l/p)-l/2 foj. g > iff 

for all t > 0. 

A similar calculation gives 

r ||a- Ve(*-")^(M|M|^-i)(r)|Ldr 

Jt/2 

(3.20) < /'^(t-r)-i/^||<,r)||^^rfr 

< C [\t- r)-i/2(l + ^)-(n/2){.-i/p)-/35/2 

Jt/2 

< ^^-(n/2)(g-l/p)-(/3g-l)/2 



for all t > 0. 
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Finally, to obtain ( ^.18 ), combine ( p. 19 ) and (|3.20|) (note that 

-(n/2(g - 1/p) - {(3q - l)/2 < -(n/2)(l - 1/p) - 1/2 

for g > (n + 2)/(n + /?)). □ 

Proof of Corollary |2.2| . As pointed out in Section 2, it follows from Corol- 
lary ^]l| that 



(3.21) tW2)(l-lM+/3/2||^(.^^)_gtA^^^.)||^^Q ^^^^ 

for each q > 1 + l/{n + (3) and every p G [1, oo]. Thus ( |3.21| ) combined with 



Propositions |2.1| and p.l| yields Corollary |2.2| . □ 
A few remarks are in order. 

Remark 3.4. If the nonlinear term in ( [L.l| ) has the form V ■ f{u) and the 
function / is sufficiently regular at zero, it is possible to improve the conclusion 
of Corollary O to 



for the multi-index 7 depending on the regularity of /. □ 



Remark 3.5. Consider Theorem p.l| and Corollary |2.2| in the context of the 
viscous Burgers equation 

(3.22) ut-u^^ + {uy 2)^ = 0, u{x,0) = uo{x), x e R. 



This is problem ( |l.l| )-( pr2| ) with n = 1, g = 2, and a = 1/2. It is well-known 
(cf. e.g. [jl9|, ^, |13|, |, |11|, [1^) that the large time behavior of solutions to this 
equation supplemented with the integrable initial condition is described by so- 
called nonlinear diffusion waves (cf. ( |1.4|) , above). If, however, it is assumed 
that Uq satisfies the conditions from Proposition p.l| with some < (3 < 1, and 
if, moreover, sup^^g ll^'^wolli is sufficiently small, the asymptotics for large 
t of solutions to the Burgers equation is given by the self-similar solutions 
AD^G{x, t) to the heat equation. 

For completeness of the exposition, we analyze problem (|3.22|) in more detail. 
Using the Hopf-Cole transformation one obtains the solution of (|3.22| ) of the 
following form 

(e*^w;o)^(x) 

3.23 u{x,t) = - \ ; 
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where as usual wo{x) = exp J^^uo{y) dyj. Supposing that uq G L^{]R), 
jjnUQ^x) dx = 0, and Uq satisfies the Miyakawa moment condition 



(3.24) 



|a;|^|Mo(a;)|(ix < oo, 



it is easy to show directly from the explicit formula ( p.23| ) that the L^-norm of 
u{-,t) decays at the rate Indeed, first note that the denominator is the 

solution to the heat equation with the datum wq and is bounded from below 
by exp(— 11^0 111). Thus, it is only necessary to bound the numerator 



G{x - y,t){wo)yiy) dy 



M 



G{x - y,t)uQ{y)exp 



uo{z) dz I dy. 



Obviously, ("^0)^ G L^{IR) and //r |l/|^|(wo)y(y)| dy < 00, since Uq has these 
properties and exp ^— Uf){z) dzj is a bounded function. Let us skip an 
easy proof that Ji[i{wQ)y{y) dy = 0. Consequently, (u^o)y satisfies the Miyakawa 
conditions, so the L^-norms of solutions to ( p.22| ) decay with the rate t~^/^. 
Finally, repeating the calculations from the proof of Corollary |3.1| yields that 
the large time behavior of solutions to ( p.22| ) is described by AD^G{-,t). Note 
that here non-smallness assumptions on uq have been imposed unlike it was 



done in Theorem 2A in the case l + l/(n + /?)<g<l + l/n. This example 
suggests that such an assumption in Theorem |2.1| is not necessary, however, 
the proof of a stronger version requires new ideas. □ 

Remark 3.6. In this paper, we limit ourselves to the case /? G (0, 1). We 
expect a completely different large time behavior of solutions to ( |l.l| )-( pr^ ) in 
/3 > 1 for the following reason. Suppose that 



(3.25) 



Uq G L}{1R^ , (1 + |a;|) dx) and / Uq{x) dx = 0. 



It is proved in [T^ that ||e*^uo||i < Ct "'^^^||Mo||Li(iR",|a;| dx) for all t > and a 
constant C; moreover. 



t 



1/2 



xuq{x) dx • VG{x, t) 



as t — > 00. 



Now, using the second order asymptotic expansion by Zuazua ||3^ (cf. also 
[0] for analogous results with more general diffusion operators and less regular 
initial conditions) of solutions to (|l.l|)-(pr^) with g > 1 + 2/n, we obtain that 
the quantity 



^1/2 



u{-,t)—\ / xuo{x) dx — a / {u\u\'^ ^){x,t) dxdr] ■ 'VG{x,t) 

\J]R" Jo JR" J 
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tends to as t — s> oo. This asymptotic result shows that the large time be- 
havior of solutions with the initial data satisfying (|3.25|) can be classified as 
weakly nonlinear in the sense of Zuazua |3^. Here, however, the first term of 
the asymptotics comes linearly from the heat kernel, but has a nonlinear de- 
pendence on the solution through a multiplicative factor (as noted by Zuazua 
in it is an open question if this factor is different from zero). Hence, 
assuming that ||e*^Mo||i < Ct~^^'^ for some (3 > 1 one should expect asymp- 
totic expansions of solutions completely different from that in Corollary |2.2| , 
specifically of the form just described. □ 



4 Nonlinear asymptotics 

The following two lemmata give the crucial steps to yield the necessary esti- 
mates of the integral equation ( ^.11 ) . 

Lemma 4.1 Let a G iR" be a fixed constant vector. There exists a constant 
C > such that for every w G L^{IR^) we have 

(4.1) \\a ■ Ve*^w7||g-/3.o. < Ct^^-^^'^w\\i 

for all t > 0. 

Proof. Using the definition of the norm in i3f and properties of the heat 
semigroup yields 

\\a-Ve'^w\\^-0,^ = sup/Z^e^^a- Ve*^M;||i 

1 s>0 

= sup//la-Ve(*+^)^w7||i 

s>0 

< C\\w\\isups('/\t + s)-'/^ 

s>0 

for all t > 0. Now, a direct calculation shows that sup^^g ■^'^''^('^ + s)~^^'^ = 
C(/5)t(^-i)/2 with C{P) independent of t. □ 

Lemma 4.2 Assume that v G Bi^''^ . Then for each p G [1, oo] there exists a 
constant C > such that 

||e*^^7||p < Cr("/2)(i-i/p)-/3/2||^||^_^^^ 

for all t > 0. 
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Proof. Standard properties of the heat semigroup e*^ and the definition of 
the norm in B^^'^ give 

for all t > and a constant C. □ 



We are ready to prove the existence Theorem |2.2| in the critical case q*. 

Proof of Theorem |2.2| . Our reasoning is similar to that in P, ^ 
Moreover, the calculations below resemble those in the proof of Theorem |2.1| 
with l + l/{n + P) < q < 1 + 1/n, thus we shall be brief in details. Recall that 
in this section we consider 



q = q =1 

which is equivalent to 



1 



n + p 



- I 







= i i + £ 




2 


q \2 2 



2 

Equip the space X with the norm 

= max{sup |L-,.oo, sup t^^/'^^^-^/^^+^Z' 

We will show that the nonlinear operator 

(4.2) Af{u){t) = e'^uo - Ta- Ve(*-")^(u|u|'^-^)(r)dr 

Jo 

is a contraction on the box 

Bii, = {ueX : \\u{t)\\^^p,^ < R and sup t("/2)(i-i/g)+/3/2 jj^^^^ ||^ < 2e} 

1 i>0 

for sufficiently large R> and a suitably small e > 0. This will be guaranteed 
provided the following estimates can be shown to hold 

(4.3) \\Af{u){t)\\^-,3,^ < ||Moi|g-/3,oo +C£'', 

(4.4) t("/2)(i-V'/)+/3/2||Ar(M)(t)||^ < C\\uo\\^-,,o.+Ce'^, 

and 

(4.5) mum-mvmw^-.^^ 

< Ce'^-Supt("/2)(l~l/,)+/3/2||^^.^^) _^^.^^)||^ 

t>0 

(4.6) __^(^)(^)||^ 

< C£''-Supt("/2)(i-i/'?)+^/2||M(-,t) -^;(-,t)||,. 



Convection-diffusion equation 



21 



For the proof of ( [4 .31) observe that ||e*'^Uo|L-^i.oo < ||'Uo|L-/3,oo. Hence comput- 



with constants C independent of u and t. 

|ig-p,oo ^ II "-u lin- 
ing the i3]~'^'°°-norm of ( [4.2|) for u G B^ sr and applying Lemma we obtain 

\\Ar{u){t)\\^-0,^ < ||e*^Mo||g-;3,oo + \\a ■ Ve(*-")^(M|M|''-i)(r)||g-/3,oc dr 
< \\uo\L-P.^ + Ce'^ f\t - r)('3-i)/V-("/2)(<?-i)-/39/2 dr. 

1 JO 

Note now that the assumptions P G (0, 1) and q = l + + guarantee 
that the integral on the right hand side is finite for any t > 0. Moreover, since 
— 1)/2 — n(g — 1)/2 — /9g/2 + l = 0, it follows that this integral is independent 
of t. Hence, estimate ( [4 .31 ) holds true. 

The proof of ( [4.4| ) is similar. It involves Lemma |4.2| as follows 

mu){t)\\, < \\e'''uo\U+ f'\\a-Ve^'~^^''{u\ur'){r)\\,dT 



(4.7) < cr("/2)(i-i/'^)-'5/lwo|lB-^- 



In this case, the conditions on (3, q imply again that the integral on the right 
hand side is finite for every t > 0. In fact, by a change of variables, it equals 
^^-(n/2)(i-i/p)-/3/2 ^ constant C > 0. Hence (U) is proved. 
The proofs of ( [4.5|) and (^4.6| ) are completely analogous. The only difference 
consists in using elementary inequality 

(4.8) \\u\u\^~^ - vlvl^-^W^ < C\\u - v\\g [\\u\\l~^ + \\v\\l-^ 

valid for all u,v e ^^(iR"). 

Finally, it follows from ( [4.3| )-( ^l6| ) that A/" : Bji ,, 5^ ^ is a contraction for 
R > 2\\uq\\i^-i3.oo and a suitably small e > 0. Hence the sequence defined 
as uoit) = e*^Mo and Un+i(t) = A/'(u„(t)) converges to a unique (in -B/j.e) 
global-in-time solution to (|1.1|) - ([L^ ) provided Uo(t) G -B_r,£, i.e. ||uo||g-;3,oo is 
sufficiently small (cf. Lemma E]2|). □ 



The proof of Theorem |2.3| requires the following result from [^, Lemma 6.1]. 

Lemma 4.3 Let w G L^(0, 1), w > 0, and Jq w{x) dx < 1. Assume that f 
and g are two nonnegative, bounded functions such that 

(4.9) fit) < git) + [\iT)firt)dT. 
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Then limt^oo g(t) = implies limt^oo fit) = 0. □ 



The next task is to prove the stabihty Theorem p.3| . 

Proof of Theorem 2^. The subtraction of equation (|2.1|) for v from the 



analogous expression for u leads to the following identity 

u{t)-v{t) = e^^{uo-Vo) 

(4.10) ~ lo^' ^^^''"'^'^ {u\u\'''^ - v\v\'^-^) (r) dr. 

Repeating the reasoning from the proof of ( [4.4|) involving inequality ( [4.8|) gives 
\H;t)-vi;t)\\, 

(4.11) < cr^^^/m-m-m (^^^/^f/^e^'/^'^'^iuo - vo)h) 

+C r(t-r)-("/^)(i-V^)-V2|m,r)-.(,r)||, 
Jo 

xf|K,r)||r + ||<,r)||ri) dr. 



By Theorem |2.2| , the both quantities 



snp&/'^^'-'/'^^-'^/'\\ui;t)\U and supt("/2){i~i/.)+/3/2||^^_^ 

t>0 f>0 



are bounded by 2£. Hence, multiplying ( [4. Ill ) by putting 

(4.12) fit) = t(-/m~m+m w^^.^t) - v{;t) ||,, 

and changing variable r = ts, we get 

(4.13) +2Ce''-' f\l - s)-W2){i-i/.)-i/2^-(n/2){,-i)-/3,/2^^^^^ 

Since (l _ s)-(V2)(i-i/g)-i/2^-(n/2)(g-i)-/3g/2 ^ ^^1(0,1) (cf. commeuts foUowlug 
inequalities (^T?!) ), we may apply Lemma obtaining /(t) ^ as t ^ oo for 
sufficiently small e > 0. This proves ( p.l3| ) for p = q. 

Next, we prove (|2.13|) for p = 1. Computing the L^-norm of (|4.10|) and repeat- 
ing the calculations from (|4.11|) and ( [4.13| ) yield 



+ C f\l - s)-l/2s-("/2){.-l)-/3./2j(^^) 
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where /, defined in ( [4.12|) , is a bounded function satisfying limt^oo fit) = 0, 
by the first part of this proof. Hence (|2.12|) and the Lebesgue Dominated 
Convergence Theorem give 



(4.14) 



hm tf'/'\\u{;t)-vi;t)\U = 0. 



t— >oo 



The next stage of the proof deals with (|2.13|) for all p G (l,oo). The calcu- 
lations from (|3.11| ) show that ||M(-,t)||oo and ||f(-,t)||oo can be both bounded 
by C't~"/^~'^/^ for alH > and a constant C independent of t. Hence, by the 
Holder inequality and ( [4.14| ) it follows that 

\H;t)-v{;t)\\p < C\H.,t)~v{;t)\\y' 

x(||n(-,t)||L-^/^+||t;(-,t)||L"^/^ 
= as t 

where we used the following inequality 



oo, 



(4.15) 



g\g\'^~^-h\hr'\<%-h\{\gr' + \h\ 
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valid for all g,h & M and q > 1. 
Finally, the proof of ( 2.13|) for p = oo involves equation (^4. 10|) and ( p. 13 



proved already for all p G [l,cx3). Standard — L'^ estimates of the of the 
heat semigroup imply that 

t^/'+^/^\\e'^iuo-Vo)\\oo < Cr/2+/5/2(^/2)-/2||e{t/2)A(^^_^^)||^ 

= C(V2)'3/le(*/^)^(«o-t;o)||i-0 



as t — > cxD by assumption (|2.12 



To study the second term on the right hand side of ([4.10|) , the integration 
range with respect to r is decomposed into [0,t] = [0,t/2] U [t/2,t]. 
Combining inequality ( 4.15| ) with estimates of the heat semi-group and the 
Holder inequality yields 



I u\u 



g-i 



v\v 



9-1 



(4.16) 



< C{t - ^)-"/2-l/2^-/3/2-(n+/3)(<ir-l)/2 



where C is independent of t and r, and /i (r) = r^/2||-u(r) - v{t)\\i is the 
bounded function which tends to as t ^ oo by ( p. 13] ) for p = 1. 
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Moreover, choosing 1/r + 1/z = 1, similar calculations lead to 



g~l 



V\V 



9-1 



(4.17) < C{t - 7-)-("/2){l-lA)-l/2^-(W2)(l-l/r)-/3/2-(n+/3)(g-l)/2 j^^^^^ 

where /^(r) = r("/2)(i-i/r-)+/3/2 _ ^(^)||^ g^jg^ ^g^^g to as t ^ oo by 
(|2.13|) . Hence, by the change of variables r = ts, it follows from (^4.16|) that 

('t/2 



1/2 



< cr"/2-^/2 / (i-s)-"/2-i/2s-^/2-("+^)('^-i)/Vi(5t)^is. 



The integral on the right hand side is finite (recall that q = l + l/(n + f3)) 
because 

2 



> -1 for /3g (0,1). 



2 2 

This integral tends to as t ^ oo by the Lebesgue Dominated Convergence 
Theorem. 

The case of the integral //^2 ■■■ d'T involves inequality (|4.17| ) with z > 1 chosen 
such that — (?7,/2)(l — 1/z) — 1/2 > —1. The proof here is analogous as in the 
last case and as such will be omitted. This completes the proof of Theorem 

□ 



5 Balance case: self-similar solutions 

In this section, we continue our analysis on the asymptotic behaviour of solu- 
tions of (|1.1|) when q is the critical exponent q = q* = 1 + l/{n + (3). Here, 



we would like to explain how Theorem |2]^ ensures the existence of a new class 



of self-similar solutions to ( |1 . 1|) and how Theorem |2.3| shows that there is a 
large class of solutions whose asymptotic behaviour in L^(iR") corresponds to 
self-similar solutions. 

Elementary calculations show that if u{x,t) is a solution to the equation 

(5.1) Mt- Af + a- V(m|m|^/("+^)) = 0, 

then so is X"'~^^u{Xx, X^t) for every A > 0. Self-similar solutions should satisfy 
the equality u{x,t) = A"'*'^m(Ax, A^t), hence choosing A = A(t) = 1/v^ yields 
a self- similar form 

(5.2) u{x,t)=r^u(-^Y 
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where U{x) = u{x, 1), x G iR", and t > 0. Substituting u{x,t) defined in ( ^l2|) 
to equation (|5.1| ) we shows the function U = U{x) satisfies the elhptic equation 

(5.3) -AU--x-VU= "^-^U + a ■ V(t/|[/|^/("+^)) = 0. 

We beheve that one can obtain solutions to ( |5.3| ) using ideas similar to those 
developed in [Q]. In that paper, Aguirre, Escobedo and Zuazua establish a pri- 
ori estimates and existence of solutions to the system 

1 TL 

(5.4) _ A/ - -a; ■ V/ = -/ + a ■ V^{f ) = 0. 

The main difference between our case and (|5.4| ) is that their coefficient for / 
is exactly n/2 which is the first eigenvalue of L = —A/ — |a; ■ V/. 
In our paper, however, we propose a completely different construction of self- 
similar solutions, based on the Cannone method 0. Let us formulate this 
result. 



Theorem 5.1 Assume that uq E Bi^'°° is a homogeneous distribution of de- 
gree —n — p. Under the assumptions of Theorem the constructed solution 
to is self-similar; hence, of the form (W^- 



The proof of this theorem follows the standard reasoning (cf. e.g. P, Section 
3] and 0, ^) and is based on the uniqueness result from Theorem 
Let us skip other details. Here, we only mention that the fractional derivative 
of order P of the Dirac delta D^5o belongs to Bi^'°°. Indeed, this follows 
from the definitions of e*^ and 6q, since c^^D^Sq = Df^G{-,t) (cf. the proof of 
Proposition |2.1|) . Hence, the self-similar form of D^G{x,t) (see ( p.l|) ) yields 
lle^^D^^olli = t-f^/^WDf^Gi-, l)||i. Finally, note that the tempered distribution 
D^Sq is homogeneous of degree —n — (3. Consequently, Theorem ^Tl] implies 
that every solution to (|1.1|)-([TT^) corresponding to AD^5q with sufficiently 
small 1^1 is self-similar. 

Now, let t~*^""^^^/^t/A(x/-\/t) denote the self-similar solution corresponding to 
the initial datum uq = AD^6q for some A E M. In the following theorem, we 
show that Ua describes the asymptotic behavior of a large class of solutions 

to (O)-(O)- 



Theorem 5.2 Let the assumptions from Theorem |^. dj hold true. Assume that 
V is the solution of ( \1.1\ )- ^T^ constructed in Theorem \2.2{ corresponding to 



26 



G. Karch k M.E. Schonbek 



the initial data vq G Bi^'°° . Let t~^'^~^^^^'^UA{x/ y/t) be the self-similar solu- 
tion corresponding to the initial datum Uq = AD^6o for sufficiently small \A\. 
Suppose that 

(5.5) \im t^/^e'\o - AD^G{-,t)\\^ = 0. 

t—^00 

Choosing e > in Theorem sufficiently small, we have 

lim t("/2)(l-l/P)+/3/2||^(.^t) _ ^-(n+/3)/2f;^( /^)|| ^ q 
>oo 

for every p E [1, oo]. 



This theorem is a direct corollary of Theorem p.3| . Recall only that, by Propo- 
sition 2A_, the limit relation in ( |5.5|) holds true if, in particular, I^vq G L}{IR"'). 



In this case, A = Jjp^„ Ipv^i^x) dx. 

Let us compare Theorem with its counterpart proved by Escobedo and 
Zuazua in |TB|, and recalled already in Introduction, formula ( |1.4| ). When 
J uq = M and q = 1 + equation (|l.lj ) has a one-parameter family 
of self-similar solutions parameterized by M. Moreover, Um describes the 
large time asymptotics of all solutions with mass M. Note that for every 
Uq G L^{1R"'), the condition / mq = M is equivalent to 



e*^Mo - MG(-,t)||i ^ as t ^ oo. 



In our case, when M = 0, the set of self-similar solutions to (|1 . 1|) with q = 
1 + l/{n-\- (3) is more complicated, however, relation (|5.5|) (or, more generally, 
(p.l2|) ) still allows us to identify solutions to (|1.1|) - (|1.2|) with the given self- 



similar large time behavior. 



6 Conclusions 



The ideas developed in this paper can be applied to other types of equations. 
As the first example, let us look at the Navier-Stokes equations for the incom- 
pressible fluid 

Ut — Au + u ■ Vm + Vp = 0, 
divM = 0, 

m(-,0) = Uq. 



It well-known (see e.g. |^) that any integrable solenoidal smooth vector field 
Uq (i.e. divMo = 0) satisfies J uq = 0. This fact motivated Miyakawa to study 
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in |27] the L^-decay of solutions to the Navier-Stokes system endowed with 
integrable initial conditions satisfying ||e*'^Mo||i < Ct^'^/^ for some < /3 < 1, 
a constant C, and all t > 0. We believe that our methods will offer some 
improvements to the Miyakawa results. 

We also expect that it will be possible to improve asymptotic expansions of so- 
lutions to the Korteweg-de Vries-Burgers equation as well as to the Benjamin- 



Bona-Mahony-Burgers equation obtained recently in |22, 23 . Some prelimi 



nary progress in this direction was already done by M. Mei in [25 
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